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ABSTRACT 


This thesis presents a comparison of three Magnetic 
Anomaly Detection (MAD) models: a cross-correlation 
detection model, a square law detection model, and a model 
referred to as the OPTEVFOR detection model. FORTRAN and 
BASIC programs for the three detection models are included 
in this thesis. The programs yield detection probabilities 
for straight line encounters. Magnetic signal values for 
the straight line encounters are an additional output. 
Plots of lateral range curves and magnetic signal values are 
presented. A discussion of the required parameters is 
included in the thesis to facilitate the use of the 
programs. The parameters that were considered in the 
comparison of the three detection models are: magnetic 
noise, aircraft and submarine headings, submarine 
displacement, and the vertical separation between submarine 


апа aircraft. 
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ts INTRODUCTION 


This thesis presents a comparison of three Magnetic 
Anomaly Detection (MAD) models. The comparison is in terms 
of probabilities of detection that were computed using the 
models. Two of the models, the cross-correlation model and 
the square law model, have been used to model sonar 
detection [Ref. 1: pp. 343-357]. The third model, referred 
toas the OPTEVFOR model, is a slant range threshold 
detection model. The results of the comparisons are 
presented in graphical and tabular form. In addition, plots 
of magnetic signals for selected lateral ranges and noise 
levels are shown. The effects of noise, aircraft and 
submarine headings, submarine displacement, and vertical 
separation are also indicated. 

The models were implemented using the FORTRAN and BASIC 
programs! that are listed in Appendix A. For Lthose 
interested in using the programs for other investigations, 
an input parameter discussion is provided in Chapter 3. To 
use the FORTRAN program, the user specifies the input 
parameters in an input file. After execution of the 


program, an output file is generated that contains 


! The programs are based on an unpublished BASIC program 
by R.N. Forrest for an H.P.- 85 microcomputer. 





Brobabilities of detection for each of the three models. Іп 
addition, magnetic signal values and magnetic signal values 
plus random magnetic noise values for one of the encounters 
generated by the program are included in the output file. 
An IBM GRAFSTAT graphical package was used to produce the 
graphics in this thesis. 

To use the BASIC program, the user must interactively 
enter the input parameter values for each encounter. After 
r [ion of the program, an optional hardcopy printout 
supplies the input parameter values and a table of detection 
probabilities for each of the three models (see Appendix A). 
Following this, lateral range curves are displayed to the 
user for immediate observation. A typical program run 
producing 21 detection probabilities for each model requires 
approximately 10 minutes of computing time on an Atari 800 


microcomputer. 





fi.” MODEL DESCRIPTIONS 


THE CROSS-CORRELATION AND SQUARE LAW DETECTION MODELS 

The cross-correlation and square law detection models 
are described in detail by Forrest [Ref. 2: pp. 33-35]. The 
models are based on the following assumptions: the noise is 
gaussian, and the signal sample points are such that 
adjacent magnetic noise samples are independent. 

The magnetic signal values, as measured by an 
aircraft's magnetometer for the cross-correlation and 
square law detection models, are the submarine magnetic 
field values in the direction of the earth's magnetic field 
at the positions of the magnetometer. The submarine's field 
sumed to bea dipole field, and the aircraft and the 
submarine are assumed to keep constant speeds and headings 
111൩൮ an encounter. 

For the cross-correlation model, a complete prior 
knowledge of the magnetic signal is required. 
Operationally, this suggests that a signal file, which 
contains a replica of the signal for each possible encounter 
situation, would be required. The model describes a 
perfect detection system with respect to the noise model 
that is used. For the square law model, a signal replica is 


Bet required. This model might be considered to describe 





the limiting detection capability for an automatic system 
that does not use information about the shape of the 


magnetic signal. 


В. OPTEVFOR MAD DETECTION MODEL 

The OPTEVFOR model is described by Forrest [Ref. 3: pp. 
7-8]. In characterizing the submarine magnetic signal as a 
Simple dipole signal, the U.S. National Defense Research 
Committee, [Ref. 4: p. 20], reports that the magnetic signal 
of the submarine "varies as the inverse cube of the distance 
БСШ the source’. In an OPTEVFOR report [Ref. 5: p. 1, encl. 
1], the results of a regression analysis on empirical peak 
to peak signal output against slant range between submarines 
and aircraft are reported. These results also suggested 
this inverse cube relationship for the magnetic signal. 
This relationship is the basis for the OPTEVFOR detection 
model. 

The model has a deterministic mode and a stochastic 
mode, each of which involves the following parameters: the 
Submarine magnetic moment (M), an Operator Recognition 
Factor (ORF), the average peak to peak magnetic noise (N) in 
the operating area, and a slant range (R). The relationship 


between these quantities is given by: 


см |1/3 


R= |------- (eqn 2.1) 
(ORF) N : 


iz 





The value of the constant c is 0.10 for M in oersted 
centimeters?, R in meters, and N in gamma. 

In the deterministic mode, detection occurs if and only 
if the aircraft's slant range from the submarine at CPA is 
less than or equal to R. This mode yields a rectangular 
("cookie cutter") lateral range curve with the probability 
of detection equal to 1 for an encounter where the slant 
range at CPA is less than or equal to R, and 0 when it is 
greater than R. 

The stochastic mode allows a more uncertain approach to 
detection by allowing a gradual rise in probability of 
detection as the slant range at CPA decreases. In this mode 
one sets the probability of detection at R equal to 50 
percent, and the lateral range curve is given by pg = D(x); 
where it is understood that © is the standard normal 
cumulative distribution function and x is determined by the 


following equation: 


(AL) R ы (еап 2.2) 
In this equation, CPA is the magnitude of the slant range 
Gistance at CPA, and R is the calculated range from Equation 
2.1. The product (AL)R represents a standard deviation. 
The value of AL can be considered to be determined by "the 
combined uncertainty and variability in the values of M, N, 
and ORF" [Ref. 3: p. 8]. Two values of AL (.20 and .01) are 


shown in Figure 2.1. If empirical data was available, the 


119 





value oí AL could be chosen to provide a best fit to the 


observed results. 


Note, 


as AL approaches 0, 


mode approximates the deterministic mode. 
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1 1 11111 


The input parameters for the FORTRAN program are all 
contained in one input file. This allows parameter values 
to be easily changed without recompiling the main program or 
subroutines. Also, with a few changes, this program could 
be altered to operate in conjunction with a larger program 
enri eld a probability of detection on an individual MAD 
TUN, 

The input parameters are divided into four areas for 
discussion. They are: (1) sample interval, (2) earth 
magnetic field, (3) submarine moments, and (4) other 


inputs. 


A. SAMPLE INTERVAL 

The choice of a sample interval is discussed by Forrest 
ВЕ. 2: pp. 27-30]. In the program, the total observation 
time in seconds over which the samples are taken is entered 
in T7. This time should be long enough to encompass a 
Ncomplete signal" at the maximum expected detection slant 
range. 

As the slant range from the submarine to the 
magnetometer increases, the distance over which a 
Significant magnetic signal is present at the magnetometer 


also increases. Figure 3.1 graphically shows the difference 


15 








in the amount of time that a signal is present for slant 
ranges of 200 meters and 805 meters. In this thesis, the 
total time for a straight line encounter is assumed to be 20 
seconds. As can be be seen from Figure 3.1, a 20 second 
interval adequately covers the significant portion of the 


magnetic signal for an 805 mmeter slant range at CPA. 


200 METERS CPA SLANT RANGE 805 METERS CPA SLANT RANGE 








i 


0.12 
| 
| 
| 
| 
| 


: ПЕ. | 
—— ` A e 4 ——À. 


.- * പി ЕЕ :. 
: } 1 T 
4 








0.08 
| 
| 
| 
{ 
1 
| 


ome dr den eases js ---4------%-----4-----ь 


: 1 
. e 


GAMMA 
GAMMA 
0.04 


0 2 
PE ക തത nd 


: : 1 j 
t ; t , 
ഇട്ട %------4------- 


. . 
— осо оо ера» о ч 0 ар олн ത്ത 9 рет отео 8 
. 


i 
i 4 "A Ape e um um faee memi- -- 
* » i 


۰ . 
دف ممم مهه جه د مک ...+ о 9 в‏ 8 سه qu qae o‏ 0 وه هوه v.‏ °" 
‹ 1 + 


SECONDS SECONDS 


Figure 3.1. Magnetic Signals for Slant Ranges of 200 Meters 
and 805 Meters. 


The time between samples is set equal to the reciprocal 
of twice the upper bandpass filter frequency of the MAD 
sensor. A value of 0.9 Hz was suggested for use by Texas 
Instruments [Ref. 6: p. 112] as an upper bandpass filter 
limit in a discussion on the effects of noise on a MAD 
system. This value yields a time interval between samples 
90755 seconds. 

The sample interval length and the false alarm rate (the 


expected number of false alarms per hour) determine the 


16 





false alarm probability. The false alarm rate (F2) is 
assigned a value of 3 based on a report by OPTEVFOR [Ref. 5: 


oa 2.1]. 


B. EARTH MAGNETIC FIELD 

Input values for the earth magnetic field intensity and 
inclination, or dip angle, may be taken from two Defense 
Mapping Agency Hydrographic Center charts, [Refs. 7 and 8 
respectively], or approximated by using a program. If chart 
values are entered, the earth field intensity must be in 
units of gamma and the inclination in decimal degrees. The 
program used to determine the intensity of the earth field 
and inclination is based on a simple dipole field model that 
is described by Forrest [Ref. 9: pp. 39-43]. 

Table III-1 displays the program output values of 
inclination in decimal degrees and earth magnetic field in 
gamma for selected geographic locations. Ira ac dug ona 
corresponding values obtained from the Defense Mapping 
Agency Hydrographic Center Charts Number 30 and Number 39 
are also displayed. The last three columns are the average 
slant range in meters at which a 50 percent probability of 
detection is obtained for the three program detection 
models. The program input parameters for these slant ranges 
were the same as the base case, except for the following 
differences: a sample interval time of 40 seconds, aircraft 


and submarine headings of 0 degrees, and a submarine 
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displacement of 7,000 tons. The correlation between the 
slant ranges, comparing the chart values and model values, 
was found to be 95 to 96 percent for the three models. This 
suggests that, even though differences exist between the 
chart values and model values, there is a high degree of 
correlation in the final output. 

A limitation to the simple dipole field model is that 
it does not give an angle of declination (variation) with 


sufficient accuracy.“ 


As a result, all headings entered 
1nto this program must be in magnetic degrees. The Phoenix 
Corporation [Ref. 10: pp. 24-25] reports on geomagnetic 
field models that can represent the earth field "with 
overall accuracies better than approximately 150-200 gammas 
in magnitude and .2 in direction of the field." This 
degree of accuracy is not needed for this program, but a 
simplified version of one of these models that provided 
satisfactory angles of declination would be beneficial if 


the program were to be incorporated into a larger model that 


utilized true headings as inputs. 


e. SUBMARINE MAGNETIC DIPOLE MOMENT 
If a submarine's magnetic dipole moment is known for the 


geographical location and the submarine's magnetic heading, 


“Private communication Irom R.N. Forrest, who 
investigated the use of the simple dipole model for this 
purpose. 
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the following values may be entered in the program: (1)P, 
its magnitude in oersted centimeters cubed, (2) A, its 
direction in decimal degrees relative to magnetic north, and 
(3) B, its depression angle from the horizontal in decimal 
degrees. If it is not known, these values must be 
calculated for a specific location and magnetic heading. A 
program is included in the main program that can be used to 
calculate these values. The program is based on a model 
described by Forrest [Ref. 9: pp. 35-38]. The input to the 
program is submarine displacement in tons. The program also 
contains coefficients which relate displacement to magnetic 
moment. The values used in the program are based on values 
cited by Texas Instruments [Ref. 6: p. 4]. 

The past history of the submarine is represented by tne 
permanent longitudinal, transverse, and vertical moments of 
the submarine (M4, M5 and M6 in oersted centimeters cubed). 
For the examples in this thesis, it was assumed that 
effective deperming had been performed and program default 


values of zero were used. 


D. OTHER PARAMETERS 
1. Headings and Speeds 
Since the simple dipole earth field model used by 
the program does not produce accurate angles of declination, 


magnetic headings are required. In addition, the headings 
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must be in decimal degrees. The input parameters for 
Submarine speed and aircraft speed are entered in knots. 
2. Noise 

The magnetic noise is assumed to be such that 
adjacent magnetic noise samples are independent. This 
assumption is based in part on the filtering that is 
performed on the magnetic signal by the processing system in 
a MAD detection sensor. The standard deviation of the noise 
in units of gamma is the value entered into 51. This value 
can be approximated from operational data by taking from 
one-fourth to one-sixth of the measured peak to peak 
magnetic n. [Ref. 2: pp. 28-29] 

The OPTEVFOR detection model incorporates a value of 
average peak to peak magnetic noise (N) in the inverse cube 
ви Calculation. In the program, the value of N is 
determined by multiplying the S1 entry by four. 

be Operator Recognition Factor (ORF) 

The ORF is the value of the ratio of magnetic signal 
to magnetic noise for which the average operator would 
detect a signal 503% of the time in the presence of 
background noise for a false alarm rate of 3 per hour. An 
ORF value of 3 was suggested for use by OPTEVFOR [Ref. 5: p. 
ШЕТ? |. 

pL NEDIStoanec Parameters 
Two parameters, R8 and N7, are used to define the 


points plotted on the lateral range curves. R8 is the 
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maximum positive value of the lateral range in meters for 
which a lateral range curve value is to be computed. N7 
represents the number of lateral range curve values that are 
to be computed from the maximum lateral range to zero 
lateral range. 

The vertical separation (2) is the sum of the submarine 


depth and aircraft altitude in meters. 
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IV: RESULTS 


Program outputs of the three models for a set of base 
case conditions are presented in this section. Outputs for 
variations from the base case are also presented. The 
lateral range of an encounter (the horizontal separation 
between the submarine and magnetometer when the magnetometer 
is at CPA) for a 50% probability of detection is used as a 
measure of comparison. Signal and signal plus "noise" 
traces for several cases are presented. The traces are 
based on the signal and noise models that are part of the 
cross-correlation and square law models. These idealized 
Signal traces appear to have the characteristics of actual 
Signal traces. This suggests that the signal and noise 


models might be used for training purposes. 


A. BASE CASE 

The base case conditions are listed in Table IV-1. The 
table is ordered in the same manner that the values are read 
into the program. An annotation of each entry is included 
for clarity. 

Figure 4-1 presents the lateral range curves for the 
base case. Points on the lateral range curves are indicated 
by the first letter of the name of the model from which they 


were derived. The slight asymmetry of the cross-correlation 


23 





detection model and square law detection model curves is 
reflective of the shape of the signals that are 'received' 
in these models. 


Table IV-1. Input Parameters for the Base Case 


6 twice the upper bandpass limit in seconds 

20-0 sampling time interval in seconds 

Ee) false alarms per hour 

0 Enter inclination il = yes; 0 = no)? 

20.0 area of operation latitude in decimal 
degrees 

50.0 area of operation longitude in decimal 
degrees 

45.0 submarine magnetic heading in decimal 
degrees 

[0-0 submarine speed in knots 

ЕШ5. 0 aircraft magnetic heading in decimal 
degrees 

220-0 aircraft speed in knots 

0 Enter submarine moment (1 = yes, 0 = no)? 

0 Enter earth field (1 = yes, 0 = no)? 

0 Enter submarine perm moments (1= yes, 
OU ono) 

4000.0 submarine displacement 

200.0 vertical separation in meters 

0.1 noise (standard deviation) in gamma 

500.0 maximum lateral range in meters 

50.0 divisions of lateral range 

EL 0 ORF (Operator Recognition Factor) 

"n2 variability factor for OPTEVFOR model 

0025 lateral range iteration number for the 


magnetic signal and signal plus noise in 

the output file 
Table IV-2 lists lateral detection ranges and 
corresponding Slant detection ranges at CPA for a 
probability of detection equal to 50 percent for the cross- 
correlation and square law detection models. An equivalent 
ORF value for each model is also listed. Due to the 
asymmetry of the lateral range curves for the cross- 


correlation and square law models, the average of the two 50 
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Figure 4-1. Lateral Range Curves of the Cross-Correlation 
(С), Square Law (S), and OPTEVFOR (O) Models for the Base 
Case. 


Table IV-2. The Lateral Detection Ranges, Slant Detection 
Ranges, and ORF's of the Three Models for the Base Case. 


Lateral Detection Slant Detection  ORF 


Range (meters) Range (meters) 
Cross-Correlation 585 907 22 
Square Law 685 714 .44 
OPTEVFOR 3158 376 3 


percent detection ranges was used as the lateral detection 


range. The equivalent ORF values for the cross-correlation 


25 











and square law detection models were calculated using the 
slant detection range values with the following equation, 


which was obtained from Equation 2.1: 


R? N 2 eqn 4.1 
For the base case, the magnitude of the submarine field (M) 


at the submarine is 6.35 x 108 orested- cm?, the noise (N) is 


.4 gamma, and the value of the constant (c) is .1. This 
suggests that, in order to detect a magnetic signal 50 
percent of the time with a false alarm rate of 3 per hour, 
the magnetic signal to magnetic noise ratio should be .21 
for an ideal cross-correlation detector and .44 for an ideal 
square law detector. 

Using the ORF values, the cross-correlation and square 
law detection models can be used to describe the performance 
of an operator. To do this, a modified value of the 
standard deviation (oc) of the input noise can be used. The 
modified value is equal to (ORF)(0)/.21 for the cross- 
correlation detection model and (ORF)(o-)/.44 for the square 
law detection model. With these modifications, the two 
models can be used to describe the detection capability of 
an operator with a specified ORF. An example of a lateral 
range curve with the modified noise standard deviation for 
an ORF of 3 is presented in Figure 4.2 for each model. 
These curves are comparable to the lateral range curve for 


the OPTEVFOR model that is shown in Figure 4.1. 
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Figure 4.2. Cross-Correlation and Square Law Lateral Range 
Curves to Describe the Performance of an Operator with an 
ORE of 3. 

The automatic MAD system manufactured by Canada's CAE 
Electronics Ltd. is expected to produce a 50 percent in- 
crease in detection slant range [Ref. 11]. Using the 
detection slant range for the OPTEVFOR model of 376 meters, 
a 50 percent improvement would yield a detection slant range 
of 564 meters. The ORF for a detection system with this 
capability would be .88. The cross-correlation and the 
square law detection models could be used to yield lateral 
range curves for a system with an ORF of .88 by using a 
EMG standard deviation equal to .88 (0)/.21 and .88 (o)/ 
.44 respectively. Figure 4.3 snows the lateral range curves 
of the two detection models with a 50 percent improvement in 
Slant range detection. Note, with the modified noise 
Standard deviations, the models are essentially equivalent 


for the cases considered. 
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Figure 4.3. The Cross-Correlation and Square Law Models to 
Describe LRC's for the CAE Automatic Detection System. 
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SECONDS SECONDS 


Figure 4.4. Magnetic Signal and Magnetic Signal Plus 
Magnetic Noise at a Lateral Range at CPA of 0 Meters for the 


Base Case. 

Figures 4.4 and 4.5 present the magnetic signal anda 
representation of magnetic signal plus magnetic noise that 
would be received under the base conditions bv a 


magnetometer with a lateral range of 0 meters and of 780 
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meters. The signal plus noise trace was generated from 
Signal plus noise values obtained by adding a signal value 
to a gaussian noise value. The gaussian noise value was 
generated by multiplying the standard deviation of the input 
noise by a pseudo normal random number from a population 
with mean 0 апа variance 1. The pseudo normal random numbers 
were generated using LLRANDOMII, a resident program at the 


Neal Postgraduate School computer [Ref. 12: p. 2.2]. 
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Figure 4.5. Magnetic Signal and Magnetic Signal Plus 
Magnetic Noise at a Lateral Range at CPA of 780 Meters for 
the Base Case. 

The magnitude of the magnetic signal shown in Figure 4.4 
is very large in comparison to the background noise. The 
peak to peak signal to noise ratio is approximately 14 to 1. 
An operator would have little difficulty identifying the 
Signal in this signal plus noise trace. 

Conversely, the magnetic signal shown in Figure 4.5 is 


small compared to the background noise. The peak to peak 
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Ба! to noise ratio is .35. The probabilities of 
detection for the lateral range of 780 meters are: .95 for 
the cross-correlation detection model, .28 for the square 
law detection model, and 0 for the OPTEVFOR detection model. 
It seems apparent that an operator would have a difficult, 
if not impossible, time in detecting this signal at a 


reasonable false alarm rate. 


ЕЖ DIEFERENT NOISE INPUTS 

The first variation on the base case shows the effect of 
different noise inputs. The standard deviation (®) of the 
peak to peak noise is the input parameter that is varied. 
Table IV-3 lists the different @ values and the 


corresponding lateral detection ranges. 


Table IV-3. The Effect of Noise on Detection Range. 


Standard Deviation Lateral Detection Range in Meters 

of Noise in 

Gamma Cross- 

Correlation Square Law OPTEVFOR 

.005 2250(2259)* 1792(1803)* 1000 (1020)* 
. 01 1832 (1843) 1446 (1460) 782 (807) 
205 1110 (1128) 868 (890) 427 (472) 
2 885 (907) 685 (714) 318 4379) 
5 SE 212205505) 382 A431) 90 (219) 


“The numbers in parentheses are the slant range 
distances in meters. The vertical separation is 200 meters. 


Figure 4.6 displays lateral range curves for the three 
models when the standard deviation of the noise is .01 


gamma. These three curves show an increase in lateral 
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detection range over the base case. Note that the asymmetry 
of the cross-correlation and square law detection models is 
more apparent in Figure 4.6 than it was in Figure 4.1. 
Figure 4.7 displays the magnetic signal (which is the 
Same as the signal in Figure 4.5) and the magnetic signal 
plus magnetic noise at a horizontal distance of 780 meters 
when the magnetometer is at CPA. The signal to noise ratio 
is 3.5. The figure suggests that a MAD operator, in this 
case, should have the ability to detect a signal at 780 


meeers lateral range with a satisfactory false alarm rate. 
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Figure 4.6. Lateral Range Curves for the Three Models with 
the Standard Deviation of the Noise Set to .01 Gamma. 
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Figure 4.7. Magnetic Signal and Magnetic Signal Plus 
Magnetic Noise with the Standard Deviation of Noise = .01 
Gamma at 780 Meters Lateral Range. 
С. DIFFERENT HEADINGS 

The headings of a submarine and an aircraft in an 
encounter have an effect on detection ranges. The effect of 
different headings was investigated using the square law 
detection model, and the results in terms of lateral 
detection ranges are presented in Table IV-4. This table 
suggests that a submarine should choose a magnetic heading 
of either East or West, and, for an encounter, an aircraft 
should also choose a magnetic heading of East or West. 

Figure 4.8 shows lateral range curves for a submarine 
heading North and an aircraft heading East. In this case, 
both the cross-correlation and square law detection model 
lateral range curves display noticeable asymmetry. The 


OPTEVFOR detection model lateral range curve is symmetric 
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Table IV-4. Square Law Lateral Detection Ranges for 
Different Submarine and Aircraft Magnetic Headings 


Aircraft Submarine Headings (magnetic) 
Headings 

(magnetic) 

0 45 90 1735 180 225 270 3185 

0 T90 650 498 640 686 640 498 650 

45 724 637 505 670 ТІН ЕЕ) 624 501 676 

90 730 654 S27 654 730 64 6 524 646 

35 712 685 >05 636 724 682 502 624 

180 085 640 498 650 700 650 498 640 

225 712 624 501 6 72 724 636 506 672 

270 730 646 519 646 730 654 509 654 

315 724 685 501 624 22403 682 505 6 36 
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Figure 4.8. Lateral Range Curves for the Submarine Heading 
North and the Aircraft Heading East. 
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but, like the curves for the other models, it shows an 
increase in detection ranges over those for the base case 
(where the submarine is heading NE and the aircraft is 
heading NW). 

The APAIR MOD 2.6 [Ref. 13: p. 83] simulation uses a MAD 
detection model that accounts for the change in a 
submarine's magnetic moment (which is dependent on changes 
in submarine heading) by using a parameter labeled DFACTR 
(degradation factor for heading). In the model, D (a 
modified slant range at CPA) determines the probability of 
detection. The value of D is determined using the following 
relation: 

D DE E DEACT TI KAY, eqn. 4.1 
where DC is the slant range at CPA and A is the acute angle 
in decimal degrees between the submarine heading and an 
East -West bearing. The probability of MAD detection is 
determined from a table of probability of detection against 
slant range. A uniform (0, 1) random number is drawn to 
determine whether or not the submarine is detected. Тһе 
average slant detection ranges (computed from Table IV-4, 
where the vertical separation is 200 meters) for submarine 
headings of North and East are 741 meters and 545 meters 
respectively. These ranges yield a value of .003 for 
DFACTR. The average slant detection range from Table IV-4 
for a submarine heading of NE is 682 meters; however, the 


slant range determined by a modified slant range of 545 
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meters and a DFACTR = .003 is 643 meters. If sin A instead 
seen TS used in Equation 4.1, then DFACTR is .265 and the 
slant detection for a submarine heading NE is 670 meters. 
Since this is only a single data point and there is no 
Sporting operational data, the modification is not 
proposed as one that should be adopted. However, this 
cursory analysis does indicate a way in which the programs 
presented in this thesis might be used by others. 

Table IV-5 lists lateral ranges for P(det) equal to 50 
percent for 3 submarine/aircraft heading combinations. The 
cross-correlation and OPTEVFOR detection model results show 
the same relationship as the results of the square law 


detection model. 


Table IV-5. Lateral Ranges for P(det) = .50 in Meters for 
the Three Detection Models. 
Submarine 45 0 90 
Aircraft 5 90 0 
Cross-Correlation 885 934 754 
Square Law 685 70 498 
OPTEVFOR 218 358 230 


For the detection ranges reported by OPTEVFOR [Ref. 5: 
p. 5.1], the effect of different headings was averaged out. 
That is, measurements were taken from the 16 possible 
combinations of the 4 cardinal submarine and aircraft 
headings in equal numbers and then averaged to yield an 
average slant detection range. But, as shown in Tables IV-4 


and IV-5, the models show significant variability in lateral 
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detection range for different submarine and/or aircraft 
headings. 

Figure 4.9 is included to show the lateral range curves 
when the submarine is headed East and the aircraft is headed 
North. These lateral range curves give the minimum lateral 
detection ranges for the different heading combinations. 
Also, for the cross-correlation and square law detection 


models, the lateral range curves are fairly symmetric. 
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Figure 4.9. Lateral Range Curves for the Submarine Headed 
East and the Aircraft Headed North. 
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D. SUBMARINE DISPLACEMENT 

The submarine magnetic dipole moment program within the 
main program is used to calculate a submarine's induced 
magnetic moments. The program is based on a model described 
Meee Oorrest [Ref. 9: pp. 35-38]. The model requires 
submarine displacement as an input. Table IV-6 displays 
results when the submarine displacement is doubled in each 
succeeding entry. 
Table IV-6.Slant Detection Ranges in Meters for 

Different Submarine Tonnages. 

Displacement Signal Slant Detection Ranges in Meters 


in tons Magnitude Cross- Square OPTEVFOR 
in oergted Correlation Law 
cm 


1000 1.59х108 590 463 DE 
2000 3.17х1 08 732 575 297 
4000 6.35х108 907 714 376 
8000 1.27х103 27 885 472 
16000 2.54x10? 1402 1099 597 
82000 . 5.08x10? 1724 1868 753 


As can be seen from column two in Table IV-6, the dipole 
moment is proportional to the displacement. Since the three 
detection models give a slant detection range that is 
proportional to the cube root of the dipole moment, doubling 
the submarine displacement should multiply the slant 
detection range by 21/3 ГБ. та: 1S confirmed Dy 
comparing the slant detection ranges between the entries in 
Table IV-6. Doubling the displacement multiplies the slant 


detection range by 1.24 for the cross-correlation and square 


D 








law detection models and, as expected, by 1.26 for the 
OPTEVFOR detection model. 

Table IV-7 lists the displacement in tons of selected 
Soviet submarines. The values were taken from Combat Fleets 
Еве World 1982/1983 [Ref. 14: pp. 602-614]. This table 


Table IV-7. Selected Soviet Submarine Displacements. 


Class Displacement in Tons 
Typhoon 25-30,000 
Delta III 10,500 3.2.50 
Yankee 8,000-9,600 
Reno. .rI 5,000-6,000 
1 I 4,300-5,100 
Charlie I 4,000-4,900 
Tango 3,000-3,700 
Boxtrot 1,950-2,400 
Whiskey 1,080-1,450 


is presented solely for the purpose of the information it 
contains. The submarine magnetic dipole moment program 
Should not be expected to give accurate estimates of these 
submarine's induced magnetic moments, since the program uses 
a value that relates displacement to magnetic moment that is 


based on submarines of World War II. 


E. VERTICAL SEPARATION 

Figure 4.10 shows three lateral range curves for a 
vertical separation of 500 meters. The OPTEVFOR detection 
model lateral range curve shows only a slight detection 
probability even when the aircraft passes directly over the 
submarine. The cross-correlation and square law detection 


model lateral range curves show an increase in lateral 
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detection range over the base case. The dip in the lateral 
range curves, for each of these models, suggests the complex 


variation of the magnetic signal with lateral range. 
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Figure 4.10. Lateral Range Curves for a Vertical Separation 
of 500 Meters. 

Table IV-8 lists the lateral detection ranges for 
different vertical separations. It should be kept in mind 
that these values are for a single geographic location; 
consequently, they may not be representative of other 


locations. Note that both the cross-correlation and square 
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law detection models lateral detection ranges increase with 


an increase in vertical separation until about 500 meters. 


Table IV-8. Lateral Detection Ranges for Different Vertical 
Separations. 


Vertical Separation Lateral Detection Range in Meters 

in meters Cross- Square Law OPTEVFOR 
Correlation 

100 804 614 360 
200 665 685 318 
300 942 720 225 
400 974 724 -- 
500 980 699 
600 974 629 
700 936 262 


“No longer attains a probability of detection equal to 50 
percent. 


A factor related to vertical separation is the effect of 
ocean wave noise on a MAD system. As the altitude of a 
magnetometer is decreased, the magnitude of the ocean wave 
noise increases. Because of the rate of this increase, for 
a given submarine and submarine depth there is a minimum 
altitude at which an aircraft should prosecute a submarine 
using MAD. Further investigation using an ocean wave noise 


model might be valuable. 
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Мы CONCLUSIONS 


This thesis has presented a comparison of three MAD 
detection models. The cross-correlation detection model, 
which models an optimum detector under the conditions of the 
detection model, yields the maximum detection range for a 
set of given conditions. The square law detection model 
does not describe an optimum detector under the conditions 
of the model and yields shorter detection ranges. In the 
stochastic mode, with an appropriate choice for the 
parameter AL that determines the standard deviation, the 
lateral range curves for the OPTEVFOR detection model become 
Similar to the other two detection models. Detection ranges 
for the OPTEVFOR detection model depend on the choice for 
the Operator Recognition Factor (ORF). With a value of 3 
for the ORF, it yields the shortest detection ranges. 
Adjusting the magnetic noise level by an amount proportional 
to the effective ORF, the cross-correlation and square law 
models can be used to describe the performance of an 
operator or an automatic detection system. 

The magnetic signal and magnetic signal plus noise 
traces appear to have the characteristics of actual signal 
traces. This suggests that the signal and noise models, 
which are the basis for the cross-correlation and square law 


detection models, might be useful for training purposes. 
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Variations ona set of base case parameters were used to 
show relative changes in the detection models. The 
parameters included: magnetic noise, submarine and aircraft 
magnetic headings, submarine displacement, and vertical 
separation. Significant results were the large asymmetry of 
the lateral range curves under certain conditions and the 
variation of the magnetic signal as shown by the changes in 
vertical separation. 

The FORTRAN and BASIC programs, along with an input 
parameter discussion, are included to facilitate the use of 
the three MAD detection models as they are implemented by 


the programs. 
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10 DIM G(290),D1¢18@),D2¢6190),K C109), X2C1000 
IS DIM 0361040) 

20 DEG 

zO PRINT "MAX FREQ"; 

48 INPUT Fi 

3Q PRINT :PRINT "MAX FREQ = ഥു 

69 Ti=1/Fi 

70 PRINT "INTERVAL TIME = "sSs:INPUT T7 
30 G=T7/2/T1 

39 H=INT(G) 

190 H=H+INT (C2*¢G—H) >) 

110 M=2*H+1 

122 ІР М) 290 ТНЕМ 70 

ISO PRINT SINT TIME = STZ 

140 Т7-ТіжМ 

120 PRINT “ADJ INT TIME = "“3T7 

160 PRINT "SAMPLE SIZE = "3; 

170 PRINT :PRINT "F/A RATE "s; 

120 INPUT F2 


190 PRINT :PRINT "F/A RATE = "SFA 
209 Pi=F2x*i(M-1)»T1/3609 
210 PRINT "РЕ = "spt 


223 PRINT ‘PRINT "INPUT DIP ANGLE (1=YES, 2=NO? "33 INPUT A 
222 IF A= THEN 260 
242 PRINT “DIP ANGLE PHI "::INPUT F 
230 GOTO 420 
260 DES :L:=76:L2=100 
279 PRINT :PRINT "LATITUDE "; 
260 INPUT L 
299 PRINT :PRINT “LONGITUDE "$; 
708 INPUT O 
210 PRINT :PRINT "LAT 3 "SL:PRINT "LON = “50 
329 F=SINCO-L2)*COS¢(L) :G=COS¢O-L2) *COS(L) sH=SINCL) 
330 U=G:V=H 
242 GOSUB 1309 
358 J=J-(20-L1):G=K*SINCJ) :H=k*COS CJ) 
252 U=G:V=F 
272 GOSUB 1200 
280 F-K:G-20:R2-(COS(L1)*SIN(JD) :Q2- COS (L1) «COS CJ2) 
230 USH:V=F 
400 GOSUB 1209 
410 FSATNC2*(SINCJ)/COSCI))) 
420 PRINT “PHI = "SF 
450 PRINT :PRINT “DIPOLE COURSE “s:INPUT C1:PRINT “DIPOLE SPEED “3: INPUT V1 
440 PRINT ‘PRINT "SENSOR COURSE "ss: INPUT C2:PRINT "SENSOR SPEED ";:INPUT V2 
ASO PRINT "DIPOLE COURSE = "3C1:PRINT “DIPOLE SPEED = "531 
462 PRINT "SENSOR COURSE = ":С2:РАІМТ "SENSOR SPEED = ";v2 
470 W12V2*SIN(CZ2) -V1*SIN(Ci) :WZ-V2«COS (C2) -V1*COS (C1) 
480 U-Wi:V-2W2 
430 GOSUB 1200 
502 С2-7:42-к 
519 DI=W0*T1+*4. 63/39 
29 PRINT "REL COURSE = ";CQ:PRINT "REL SPEED = ";W@ 
538 PRINT :PRINT "INPUT DIPOLE MOMENT (1=YES, 0=N0) "5: INPUT AA 
S40 IF AA= THEN 550 
550 PRINT :PRINT “MAGNITUDE P ";:INPUT P:PRINT “HOR ANGLE w "3: INPUT A 
SEQ PRINT :PRINT “VERT ANGLE OMEGA "3:INPUT 8 
579 GOTO 240 
SSQ PRINT :PRINT “INPUT ZARTH FIELD (1=YES,0=N0)> "3: INPUT AA 
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590 
eod 
612 
6-9 

292 
543 


= 
sl 


SE2 
679 
558 
630 
T = 
708 
712 
720 
759 


= 
v) 


748 
750 
768 


750 
738 
200 
2102 
az 

343 
ss 

Sea 
378 
acu 
390 
382 
319 
320 
323 
350 
249 
332 
368 
364 
ЗЕБ 
379 
374 
375 
358 
2320 


IF AA= ТНЕМ 628 

PRINT ‘PRINT "EARTH FIELD “S$: INPUT El 
GOTO 623 
E1=70800/SQR(2*COS CF) *COS(FI3 +1) 

PRINT :PRINT "EARTH FIELD = “SEL 
М4=0:М5=0:М6=0 


PRINT sPRINT “INPUT PERM MOMENTS (1=YES:; 2зМ02 “ss INPUT AA 
IF AR=0 THEN 590 


PRINT :PRINT "LONG MOMENT "3: INPUT MASPRINT “TRAN MOMENT “$s: INPUT MS 
PRINT “VERT MOMENT "5: INPUT MG 


PRINT :PRINT "LONG MOMENT = "MA: PRINT “TRAN MOMENT 3 "3:MS:PRINT 
"MG 


K1= 7,38 => S:KZ=1.6 
PRINT ‘PRINT “DISPLACEMENT "3 


“VERT MOMEN 


INPUT Ni 

PRINT ‘PRINT “DISPLACEMENT = ";NI 
№41 =41 

MISE 1 Ak TANI1ASINCF) :MI=MO+M6 


MS=E1*COS (CF) #NL*CK1*COS (C1) *COS (C1) +K2I*SINCCLI*SINCCL) 2 
MizMAnSINCCI2-MS«CO0S CC) :MZsMAmRCOSCCI)-MEXSINCCZI) 
M7zE&i*€"(COSCEOwWcKi-RHZO*NIM*SINCGCCIÓ)«*COSCCID) 

M1i=97+m1 :M2=MS+M2 

UsMi:VzMZ 

GOSUB 1229 

ടവ്‌ വേട: ഗ്മ, 

GOSUB 1200 

P=k:3B=J 

PRINT :PRINT "P = “SP2PRINT “м з "sQ:<PRINT "OMECA = "3B 
Х=Р1 

505ЧВ 1422 

Убау 

М7-Мж(1-2/2/МыУжчдҢ0НБК(2/2/М22А2 

PRINT SPRINT "VERT SEPARATION "35: INPUT Z 

PRINT :PRINT “VERT SEPARATION = "353Z 

PRINT :PRINT "NOISE "; 

INPUT S1 

52:5 1 м4 

PRINT :PRINT "NOISE -» "5S1 

PRINT : PRINT “MAX LATERAL RANGE “$: INPUT RS 

PRINT :PRINT “NUMBER OF INCREMENTS “$s: INPUT М7 

PRINT :PRINT "MAX LATERAL RANGE = “SRS: PRINT “NUMBER OF INCREMENTS = "N7 
PRINT :РВІМТ “ QRF “ss INPUT ORF 

PRINT ‘PRINT " QRF = "ЗОВЕ 

D425RG/N7:NB2Z*«N7 

PRINT :PRINT " ALPHA "::INPUT AL 

PRINT :PRINT " ALPHA = "SAL 

LI=-RE 

FOR E= TO N8 


1299 XU=L_2:X0(E)=L3 
1212 GOSUB 1622 
1020 L3=12+D4 

1828 GOSUB 1500 


71040 
1244 
10502 
1060 
1065 


NEXT E 

GOTO 123۵ 

GRAPHICS 3:COLOR 1 
XX2 INTCZ10/N8) 
xaca»-ua 


1070 FOR IsQ TU N8 


1089 
1990 


ന നാമാ 
DIC IOSINTCC1—-D2C 1) ടമാ 


1895 DSCIDsINTCC1-DSC1I)) 160) 


11002 


XQOCI-122X0CIO-XX 


1110 NEXT I 

1123 PLIT XØ): D19) 
115@ FOR [=1 TO NS 

1140 DRAWTO X@CI): DICI) 
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1159 
1158 
1170 
1158 
1138 
1092 
1134 
1136 
1138 
1204 
ES 

1 22 
1235 
220 
1510 
1528 
22 

155 

1240 
125 

1550 
1270 
1550 
1230 
1429 
1419 
1420 
2420 
1440 
1450 
1460 
1470 
1458 
1430 
1500 
1525 
1512 
1815 
Io 

е 
IO 
23255 
540 
1545 
1554 
1564 
1605 
1619 
1620 
1622 
1640 
1642 
1645 
1645 


NEXT I 


PLOT x0(0>,D2(0) 
FOR 1൭ TO NS 
DRAWTO XOLI):DZLI) 


NEXT I 


PLOT Х0(22.05(00 
FOR I=1 TO N8 
DRAWTO XUCI), DSCI) 


NEXT I 
PRINT 


ED FOR X FROM SAET "TRE 


GOTO 1610 


PRINT 


“FOR HARD COPY ENTER '1'"::INPUT CC 


IF CCT THEN GOTO 1052 


LPRINT 
LPRINT 
ЕРЕТМТ 
БӘНІМТ 
LPRINT 
ERRTNT 
ERRINT 
CORINT 
CORINT 
CPRINT 
LPRINT 
TIT 
ESRENT 
LPRINT 
LPRINT 
LPRINT 
ERRINT 
EPETNT 
Ee RINT 
ECS INT 
PRINT 
ERRINT 
LPRINT 
BERTNT 
LPRINT 
ЕРЕ 
LPRINT 
L2=-R8 


"MAX FREQ a "F1 

"ADJ INT TIME = ":T7 
‘SAMPL SIZE = "SM 

"F/A RATE = "3F2 

ട്‌ പാട്‌ 

"LAT =a "; 

"LON = "; 

൦7, РЕ 

‘DIPOLE COURSE = “set 
“DIPOLE SEED SVL 
"SENSOR COURSE s "C2 
SENSOR ЫР = ба уг 
5.55 ആ CO 

“REL SPEED ="3WQ 

നി) FIELD 3 “EL 

"LONG MOMENT = "sma 
"TRAN MOMENT = "3S 
"VERT MOMENT =» "МБ 
"DISPLACEMENT = “SNNL 

"р = “sp 

" vé = Е 

"СМЕСА = "ЗВ 

“VERT SEPARATION = "327 
"NOISE а "561 

"MAX LATERAL RANGE = “SRS 
"NUMBER OF INCREMENTS = “SN7 
sLPRINT “LTR RNG PD (GO) 


FOR I=@ TO N8 


EPRINT ES," 


2112" ШЕ оя 


LS=LS+D4 


NEXT I 


GOTO 1033 


PRINT 
END 


“END “ 


U=xQ:V=Z 

GOSUB 1900 

D=J:HMQ9=K 

RH=(Q. L*P/ CORF*S2) ) 20. 5575 
SIGS RH 

X=3¢RH-HO)/SIG 


GQSUB 1390 


1645 


1647 
1650 
1558 
1670 
1680 
1530 
1700 
171@ 
1729 
1720 
1740 


DS (ED =Y 


ഈ ടാ 


ടി 


“sDS¢T) 


30=C0S (B)*C0S (C3-9) :3BQ=Z0S (D) COS (B)+*SINCCI-A)-SINCDIASINCD) 


NQ==(SIN CD) COS (B)+*SIN(CO-9) ) —-COSCD) *SINCBD 


8 1 ൭05 നാ ഥുട മാ : 1 ഥുട നാ «COS CF) *SINCCQ)-SINCD)*SINCF) 
N1==(<SINCDI=COS CF) HSINCCH) I -COSCDIHSINCF) 
KL=P/10/HOAZ 

a4232*980+*81-J0*J1-N0+*N1:41=.3:*(NO+*B1+80+*N1) :A9=2+*N0Q+*N1-80+*91-J0D+J1 


50-0 


FOR IsQ 


te TL 


бз«<1-«<М-12/202ж03 


Q=3/HQ 


54 





1/20 
1762 
1770 
1750 
1730 
1500 
1810 
Sa 

1950 
1840 
1350 
1960 
1270 
1550 
1930 
1300 
1305 
1907 
1908 
1909 
1910 
7315 
1320 
1958 
1342 
1850 
1960 
1978 
1280 
1320 
000 
2210 
2022 
2025 
2050 
20702 
213 
2030 


0=1/(1+Я*ж@) ^3. 5 

ടേ ഥെ ലമാ жо 

GCI) =6:SO=S8+6+«6 

NEXT I 

RE TURN 

ചട മാ 

KCED KZ 

VE-VbE-Ki*SQROSO»/S1 

LO0=kK1+kK1*S50/ (S1:+S1):A3=M+L0:B3=1+L0/ (M+L0) 
V9=-SQR(2*V7 /B3)+SQR(2*AS/BI—-1) 

X=VE:GOSUB 1992 

DICE)mY 

X=V9:GOSUB 1992 

DICE =Y 

RETURN 

K=SORCU*U+V*V) SIF K=O THEN J=@: RETURN 

UK=U/K 2 VK=V /K 

IF UK)@. 999999 AND VK)Q0. 939999 THEN J=0:RETURN 

IF ШК) @. 3999999 THEN J=-ATNCVK/SQRC—-VK*VK+1) 9490: RETURN 
IF VK) @. 999999 THEN J=0:RETURN 

MM=ATN CUK/SQR ¢-UK*#UK+1)) S J=—ATNCVK/SORCH-VK*VK41394902IF MM(O THEN J=280-J 
RETURN 

Y=XsIF X)Q@.5 THEN Y=1-Y 

Y=SQRC(LOGCL/Y/Y)) 

G0=2. 513517:61=WY. 302833 :62=8. 210229 

H1=1. 432732 :H2=Y, 1829269:H3=1. 308E-3.3 

YzY—-CGGO-* Y» (G14G2wYo2/C1i- Y *CHic Y *CHZ-HZ*Y22) 

IF X)0.5 THEN Yz-Y 

RETURN 

YzX:IF XQ THEN Yz-Y 

Ыші/(1%2.2216412жУ9 

൫1 ൭. 313383153:2102=>-8, 336367232 :03=1. 731477393:Q4=-1. 82122559:05=1. 22027442 
IF Y>)24,.23 THEN Y=0:GOTO 2070 

PI=3. 14159265 

Y=EXP (-(Y*Y/ 23) /SORCZAP TINA (QU WR CE WR CA ZUR CHANG III) 
IF XQ THEN Y=1-Y 

Үш(ІМТ(10000жҮ22/10000 

RETURN 
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MAX FREQ = 1.3 
ADT TNT TIME = 20255955955 
SAMPLE SIZE = 37 

F/A RATE = 3 

PF = 0. J16GE6E6E6 

LAT = 20 

LON = 60 

PHI = 39. 409756973 

DIPOLE COURSE = 45 

DIOLS SPEED a TIA 

SENSOR COURSE = 315 

SENSOR SPEED = 220 

REL COURSE 3 SI. S7499 
KEL STEED a OS 
EARTH ЕТЕД ട്ടം 


LONG MOMENT = Y 
TRAN MOMENT = @ 
VERT MOMENT - 0 


DISPLACEMENT - 4000 

P = 674634292 

м ЕО! 

ОМЕСЯ = 27. 11175179 

VERT SEPARATION = 220 
NOISE = 9.1 

MAX LATERAL RANGE = 1509 
NUMBER OF INCREMENTS = 15 


LTR RNG PDCCC) PD (SL) PR CORT? 
-1500 9. 9344 2. 0133 2 
-1400 В. 2536 0.016€ Q 
- 1300 0. 9945 0.0131 Q 
-1200 Q. 1.291 0.0203 2 
-1100 022212 2. 02665 2 
-1000 9. 27692 0. 0403 g 
-200 0. 6.399 0.3801 2 
-809 0.9191 2259 Q 
-700 2. 399 0.6742 2 
-Е00 1 а. 3975 2Е-04 
-500 1 1 0.0122 

-4Q00 1 1 0. 1505 

-300 1 1 @. 567 

-200 1 1 0.5565 

-100 1 L 0.9774 

a 1 1 в. 3535 

100 1 1 0.9774 

200 1 1 0.5565 

200 1 1 Q. 567 

400 1 1 Y. 1603 

590 1 a. 3333 0.0152 
600 0. 3937 Q. 7654 2Е-д24 
780 0: 3157 0.2199 2 
590 9.5775 0.0677 2 
300 @. 3014 0.0221 2 
1900 0.1644 0.0226 2 
1100 0. 1003 0.0154 2 
1200 0.0502 മ. 2165 2 
1200 2. 2508 2. 2126 ð 
1400 2. 0405 9.0121 2 
1500 2. 254 0. 90145 Q 
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